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1.  INTRODUCTION 


CH2NNO2,  or  methylene  nitramine  (MN),  is  tiiought  to  be  one  of  the  primary  decomposition  products 
of  hexahydro-l,3^-trinitro-l,3,5-triazine  (RDX)  (Schroeder  1985a,  1985b).  MN  has  never  been  isolated 
and  is  believed  to  decompose  spontaneously  after  formation  through  concerted  channels  only  (Zhao, 
Hintsa,  and  Lee  1988).  Scission  of  the  weakest  bond  in  MN  (the  N-N  bond)  is  not  observed  (Zhao, 
Hintsa,  and  Lee  1988),  which  is  a  suiprising  result.  One  explanation  could  be  simply  that  energy  required 
for  reaction  through  concerted  channels  is  substantially  lower  than  the  scission  reaction.  Another 
explanation  is  that  there  is  anomalous  dynamical  behavior  in  the  system.  We  wish  to  perform  molecular 
dynamics  simulations  of  the  decomposition  of  MN,  but  to  do  so  requires  an  analytic  potential  energy 
surface  (PES)  that  accurately  describes  the  system  and  reaction  charmels.  The  development  of  such  a  PES 
is  the  focus  of  this  report. 

The  only  characterizations  of  MN  come  from  theoretical  predictions.  Thus,  we  are  limited  in  the 
information  that  we  can  use  to  construct  a  potential  energy  function  of  this  molecule.  We  based  the  PES 
of  MN  that  we  developed  in  this  work  on  the  ab  initio  multiconfigurational  (MC)  SC!F  and  multireference 
(MR)  Cl  electronic  structure  calculations  of  Mowrey  et  al.  (Mowrey  et  al.  1990;  Mowrey  impublished). 
In  that  work,  they  examined  what  they  believed  to  be  the  primary  decomposition  channels  for  MN: 

CH2NNO  H2CN  +  NO2  (I) 


and 


CH2NNO2  HCN  +  HONO.  (H) 

They  calculated  the  stmcture,  relative  energies,  and  frequencies  of  equilibrium  MN,  the  five-centered 
transition  state  for  HONO  elimination,  and  the  products  for  the  two  reactions  (Mowrey  et  al.  1990). 
Additionally,  they  calculated  points  along  the  reaction  path  from  the  five-centered  transition  state  of  (II) 
leading  to  reactants  (MN)  and  products  (HONO  +  HCN)  (Mowrey  et  al.  unpublished).  The  stractures  of 
the  planar  equilibrium  MN  and  five-centered  transition  state  leading  to  HONO  +  HCN  are  shown  in 
Figure  1.  After  adjusting  for  zero-point  energy  effects,  basis  sets  and  different  levels  of  theory,  they 
estimate  that  the  activation  energy  for  (II)  is  31  ±  4  kcal/mol.  They  also  estimate  the  N-N  bond 
dissociation  energy  for  (I)  to  be  35  ±  4  kcalAnol.  Although  these  estimates  suggest  that  reaction  (II)  is 
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energetically  favored,  the  steric  effects  associated  with  reaction  (H)  would  make  reaction  (I)  entropically 
favored. 

Mowrey  et  al.  (1990)  have  shown  that  the  activation  energy  for  (II)  is  lower  than,  or  at  least 
comparable  with,  the  N-N  bond  dissociation  energy.  The  analytic  PES  that  we  report  here,  based  on  the 
Mowrey  et  al.  calculations  (Mowrey  et  aL  1990;  unpublished),  can  be  used  in  molecular  dynamics 
simulations  to  investigate  whether  MN  exhibits  dynamical  behavior  that  would  preclude  reaction  (I). 


2.  POTENTIAL  ENERGY  SURFACE  (PES) 


2.1  Analytical  Form  of  the  Potential.  The  total  potential  energy  expression  is: 


2  2  K  2 


^NN  X)  ’^CH(i)  +  Y,  ^NO(i)  *  ^HCH  *  Y  '^NCH(i) 
i*l  i*l  i*l 


2 

’^CNN  Y  ’^NNO(i)  Y  '^CNNO(i)  Y  ''^NNCH(i)  *  ''^NNOj 
i*l  i*l  i*l 

*  (^ONO  ^HOn)  *  ^HONO  *  *  ^Pro 


Stt 


17  r 


i«14 


^HONO  • 


^  •  %ONO  *  “  Spro), 


(1) 


where  V®  is  the  Morse  function  that  describes  bond  stretches,  V**  is  a  harmonic  oscillator  function  that 
describes  bending  motions,  and  V“,  V',  and  Vjjqno  are  functions  of  cosines  describing  out-of-plane 
wagging  and  torsional  motions. 

These  simple  functions  (V^,  V**,  V“,  and  V',  defined  in  subheadings  in  Table  1)  are  not  flexible 
enough  to  describe  correctly  the  system  at  all  regions  of  the  PES  if  constant  parameters  are  used.  Many 
of  the  parameters  for  these  functions  must  have  significantly  different  values  at  the  product  geometry  than 
at  the  reactant  geometry.  For  example,  the  C-N  bond  in  the  product  HCN  is  much  stronger  than  the  C-N 
bond  in  the  reactant  CH2NNO2.  Therefore,  constant  parameters  for  the  C-N  interaction  terms  would  not 


3 


Table  1.  Potential  Energy  Functions  and  Parameters 


=  D{1  -  exp(-alr  -  rj)}^  -  D 
CN 


Parameter  Type 


0.847500 

-0.151844 

0.060691 


3.232720 

0.773610 

2.497460 


5.900000 

2.710000 

1.149500 


-0.666680 

0.000000 

-0.760780 


Parameter  Type 


0.759000 

1.229500 

1.062930 


4.136200 

0.975370 

2.529080 


0.000000 

0.000000 

1.441141 


Y[CH(1)]  -  Ya(Si4  +  Sjj)  +  Yg(Sig  +  S17), 
YfrHnn  =  ^15)  ■*"  ^a(^16  "*■  ^17) 


-1.321150 

-0.164100 

-0.787890 


=  D,  a,  or  r. 


Parameter 


P(+) 

P(+) 

Constant 


NO  D  =  -  DProduct) 

Y[N0(1)]  =  Ya(Si4  +  Sjg)  +  Yg(Si5  +  Sj7), 
_ '^fNO(2)1  =  +  Ya(Si5  +  Sj7) 


2.168000 

1.069000 

1.073393 

0.365500 

0.128000 

0.002680 


2.760820 

1.942640 

1.965080 

0.263500 

2.685560 


0.000000 

0.000000 

4.336000 

1.895000 

1.065370 


-0.652960 

-0.952120 

-0.970880 

0.000000 

-1.042670 


jjProduct)  S 


Parameter  Type  X 


Y  =  a,  or  r. 


Constant 

P(+) 

P(+) 

P(-) 

P(-) 

P(-) 

P(+) 


6.614000 

0.488000 

0.505000 

0.048115 

0.675500 

0.101000 

0.037044 


2.766320 

2.500000 

2.591180 

3.101050 

2.500000 

2.648640 


3.451000 

1.350000 

1.303440 

4.427000 

2.360000 

1.166310 


-1.229720 

0.000000 

-0.625800 

-0.761420 

0.000000 

-0.978410 


Y  —  jy^eactant  ^Reactant 


Parameter  Type 


Constant 

Constant 

P(-) 

P(+) 

Constant 


^OH-NN 

*^H-NN 


4.624000 

2.293640 

2.132500 

2.179910 

0.954106 


2.604230  0.000000  -1.252450 
3.127710  2.296000  0.000000 


Table  1.  Potential  Energy  Functions  and  Parameters  (continued) 


V®  =  k  (0  - 
HCH 

©e)" 

Parameter 

Type 

X 

A 

B 

C 

Ro 

k 

P(+) 

^H-NN 

0.641699 

2.851970 

0.000000 

-0.000020 

©e 

Constant 

119.91592 

CNN 

Parameter 

Type 

X 

A 

B 

c 

Ro 

k 

P(+) 

*^H-NN 

1.712481 

2.797200 

0.000000 

-0.839420 

©e 

Constant 

115.72031 

ONO 

Parameter 

Type 

X 

A 

B 

c 

Ro 

k 

P(-) 

^H-NN 

1.337147 

4.363370 

3.530805 

0.000000 

©e 

P(+) 

^H-NN 

6.539728 

3.084920 

111.38000 

-0.752150 

HCN  YfNc 

H(l)]  =  Ya(^14 

+  Sis)  +  Yb(S 

16  +  ^1?)’ 

Y  =  k  or  ©e 

^rNCH(2)1  ~ 

+  Sis)  +  Ya(S 

16  +  ^1?) 

Parameter 

Type 

X 

A 

B 

c 

Ro 

P(+) 

^H-NN 

1.198696 

3.313430 

0.000000 

-0.743490 

©eA 

Constant 

124.12487 

P(+) 

^H-NN 

0.799410 

3.882660 

0.000000 

-0.373060 

®eB 

P(-) 

*^H-NN 

32.020394 

2.128550 

115.95921 

-1.029180 

NNO  YfN, 

'40(1)1 

+  Sig)  +  Yg(Sis  +  S17), 

Y  =  k  or  0g 

1 

+  Sig)  +  Yy^(Sis  +  S17) 

Parameter 

Type 

X 

A 

B 

C 

Ro 

kA 

P(+) 

*^H-NN 

1.814668 

1.765720 

0.000000 

-0.888670 

©eA 

Constant 

121.07576 

^B 

P(+) 

*^H-NN 

1.264759 

2.549750 

0.000000 

-0.479980  ' 

©eB 

Constant 

114.46478 

HON 

Parameter 

Type 

X 

A 

B 

c 

Ro 

k 

P(-) 

^H-NN 

1.164030 

2.676090 

0.000000 

-1.002470 

©e 

Constant 

107.89000 
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Table  1.  Potential  Energy  Functions  and  Parameters  (continued) 
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'^HONO  ®  ajCosCjx) 

j-O 

Parameter 

Type 

X 

A 

B 

c 

Ro 

Constant 

0.210699 

h 

Constant 

0.024922 

32 

Constant 

-0.192175 

as 

Constant 

-0.011035 

a4 

Constant 

-0.005143 

^5 

Constant 

-0.000555 

V(t)  =  k[l  - 

Cos2(x)] 

Parameter 

Type 

X 

A 

B 

c 

Ro 

k(CNNO) 

Constant 

0.067100 

k(HCNN) 

Constant 

0.980000 

V(©)  =  kCos2(©) 

Parameter 

Type 

X 

A 

B 

c 

kCNNOj) 

Constant 

1.260581 

kCHjCN) 

Constant 

0.535061 

Attenuation  Factors 

Parameter 

Type 

X 

A 

B 

c 

Ro 

^NCKl) 

P(-) 

*^NO(l) 

0.500000 

5.000000 

0.000000 

3.500000 

SnO(2) 

P(-) 

^NO(2) 

0.500000 

5.000000 

0.000000 

3.500000 

Sjj 

P(-) 

Rnn 

0.500000 

20.00000 

0.000000 

2.500000 

Spro 

P(+) 

*k)H-NN 

0.500000 

20.00000 

0.000000 

-1.370360 

^OH-NN 

P(+) 

25.00000 

6.000000 

0.000000 

2.350000 

correctly  describe  this  interaction  at  all  regions  of  the  PES.  Also,  some  interactions  do  not  exist  at  all 
regions  of  the  PES  and  should  be  attenuated.  Examples  of  these  are  the  CNNO  and  HQW  torsional 
motions,  as  well  as  the  CNN  and  NNO  bends  and  NNO2  out-of-plane  wag.  We  have  incorporated  the 
required  flexibility  into  equation  (1)  by  using  switching  functions  either  to  attenuate  terms  or  to  vary 
parameter  values  appropriately  at  the  different  regions  of  the  PES.  The  forms  of  most  of  the  SNvitching 
functions  and  the  parameters  are  either 
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P(+)  =  A{1  +  tanh(B[X  -  Rq])}  +  C 


(2) 


or 


P(-)  =  A{1  -  tanh(B[X  -  Ro])}  +  C  .  (3) 

where  X  denotes  the  appropriate  independent  geometric  variable.  P(+)  and  P(-)  vary  continuously  in  the 
range  of  (2A  +  C)  to  C.  The  values  of  A,  B,  and  C,  as  well  as  the  definition  of  the  independent  geometric 
variable  for  the  switching  functions  and  parameters,  are  given  in  Table  1. 

One  independent  geometric  variable  used  to  vary  several  of  the  parameters  in  equation  (1)  is  defined 
as: 


^H-NN  =  • 


(4) 


Rj^jsj  is  the  N-N  bond  distance,  and  Rq^  reflects  a  weighted  contribution  from  each  of  the  four  possible 
0-H  pairs  in  the  molecule.  The  pairs  are  denoted  as  q,  i  =  14-17  and  correspond  to  O-H  pairs 
0(1)-H(1),  0(2)-H(l),  0(1)-H(2),  and  0(2>-H(2),  respectively.  Rqh  is  defined  as: 

Roh  =  X) 

i»14 


where 
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The  switching  fimction  Soj^-nn  was  added  to  equation  (4)  to  eliminate  the  (Roh“^nn)  dependence  at 
regions  of  the  PES  where  both  Rqjj  and  Rj^j^  are  very  large. 

The  description  of  Rqjj  ensures  that  unless  two  or  more  H-0  distances  are  equal,  the  Rqh  description 
is  defined  by  the  shortest  H-O  bond.  If  any  0-H  distances  are  equal,  the  contribution  to  Rqj^  is  divided 
evenly  among  the  equal  H-0  bonds.  Rqh  is  also  used  as  the  independent  variable  in  the  Morse  function 
for  the  H-0  stretching  interaction. 

The  term  to  describe  HON  bending  motions  consists  of  weighted  contributions  from  all  possible  HON 
angles  (denoted  0;,  i  =  14-17,  where  the  H-0  bond  in  the  HON  angle  corresponds  to  rj  as  defined 
previously). 


V 


b 

HON 


17 


VE 

i»14 


.HON 

1 


(7) 


where 


-  S,4  (1  *  2S,J) 
S™"  =  S,5  (1  ♦  2S,4) 
816°"  =  Sjs  (1  *  2S„) 
8,",°'"  »  S„  (1  +  2Si6)  . 


The  model  PES  allows  hydrogen  migration  in  HONO;  features  of  the  saddle  point  for  this  hydrogen 
migration  are  discussed  in  the  following  section.  This  fonn  of  the  HON  bending  potential  describes  the 
contribution  of  the  two  possible  HON  bends  in  HONO,  which  have  an  equal  contribution  when  HONO 
is  in  the  saddle  point  configuration  for  the  hydrogen  migration  reaction. 

The  attenuation  function,  Shqno*  defined  as 

^HONO  "  ^NO(l)  (^14  ^15)  %0(2)  (^15  ^1?)  •  (9) 
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attenuates  terms  upon  the  secondary  decomposition  of  HONO.  This  function  also  modifies  the  H-0 
interaction  to  correspond  to  that  of  the  diatomic  upon  HONO  decomposition. 

The  functional  fom  of  the  parameters  for  some  of  the  stretching  and  bending  interactions  are  a  bit 
more  complex  (see  Table  1)  than  just  P(+)  or  P(-)  as  defined  in  equations  (2)  and  (3),  for  the  following 
reason:  The  five-centered  transition  state  leading  to  formation  of  HONO  +  HQ4  involves  formation  of 
an  H-O  bond  while  breaking  a  C-H  bond.  Because  the  CH2NNO2  molecule  has  two  oxygens  and  two 
hydrogens,  there  are  four  possible  H— O  pairs  that  could  form  the  HO  bond  in  HONO  (as  one  of  two  C— H 
bonds  is  broken).  The  potential  energy  function  must  be  unbiased  and  invariant  to  any  of  these  bond 
formations  and  scissions;  i.e.,  HONO  +  HCN  must  be  formed  in  the  same  manner  regardless  of  which 
H-0  pair  forms  the  bond.  Additionally,  the  N-0  interaction  in  HONO  is  dramatically  different  depending 
on  whether  the  oxygen  in  the  NO  moiety  is  terminal,  bound  to  hydrogen,  or  is  a  free  diatomic.  The 
function  must  be  able  to  discriminate  in  an  unbiased  fashion  among  all  possible  bonding  situations.  The 
added  complexity  of  these  functions,  as  well  as  inclusion  of  the  various  switching  functions,  accommodate 
these  situations.  Any  modifications  of  parameter  functional  forms  are  given  in  Table  1. 

2.2  Features  of  the  PES.  The  chosen  parameters  ensure  that  the  model  reproduces  measured 
properties  of  the  products  and  reaction  endothermicities  as  shown  in  Figure  2.  The  classical  barrier 
heights  for  (I)  and  (II)  differ  by  only  2  kcalAnol,  with  (I)  being  lower  in  energy  at  35.0  kcal/mol. 
Geometric  parameters,  harmonic  vibrational  frequencies,  and  relative  energies  of  the  reactant,  the  products 
and  the  transition  state  are  given  in  Table  2,  along  with  the  experimentally  measured  or  theoretically 
determined  information.  The  relative  energies  cited  in  Mowrey  et  al.  (1990)  are  not  included  in  Table  2, 
because  the  cited  values  are  estimated  activation  energies,  and  include  corrections  for  various  levels  of 
theory,  different  basis  sets,  and  zero-point  energy  effects.  The  only  uncorrected  ab  initio  energy  value  to 
which  we  rigidly  adhered  in  fitting  the  PES  was  the  MRCI(DZ)  energy  difference  between  CH2NNO2  and 
HONO  +  HCN  (-27.9  kcal/mol)  in  order  to  establish  the  relative  energy  position  of  the  reactant. 

The  only  information  available  about  the  reactant  and  saddle  point  leading  to  formation  of 
HONO  +  HCN  are  the  ab  initio  data  of  Mowrey  et  al.  (1990;  unpublished),  which  includes  stmctures  and 
matrices  of  energy  second  derivatives  for  the  two  stationary  points,  as  well  as  several  structures  and 
energies  along  the  reaction  path  leading  from  this  saddle  point  The  most  complete  set  of  information 
provided  by  Mowrey  et  al.  (1990;  unpublished)  are  CASSCF/DZ  structures  and  second  derivatives  which 
we  have  used  for  fitting  the  PES.  The  stmctural  parameters  of  the  PES  at  equilibrium  are  in  exact 
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CHoN  +  NO 
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Table  2.  Structural  Parameters,  Zero  Point  Energies,  and  Relative  Energies' 
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Distances  given  in  A,  angles  given  in  degrees,  and  energies  given  in  kcal/mol.  ^  Chase  et  al.  (1985). 

^  Herzbcrg  (1945).  ^  Huber  and  Herzberg  (1979). 


agreement  with  the  ab  initio  values.  The  transition  state  structure  predicted  by  the  model  PES  is  in 
extremely  close  agreement  with  the  ab  initio  structure,  with  the  largest  disagreement  in  a  structural 
parameter  being  less  than  0.5%.  The  structural  parameters  predicted  by  the  model  PES  for  the  products 
(H2CN  +  NO2)  and  (HO  +  NO  +  HCN)  are  in  fair  agreement  with  the  experimental  and  ab  initio  values; 
however,  no  effort  was  made  to  incorporate  into  the  model  PES  the  small  structural  changes  as  these 
products  are  formed.  Differences  in  geometries  between  our  model  PES  and  the  experimental  and 
ab  initio  information  for  HONO  +  HCN  are  due  to  the  use  in  the  model  PES  of  previously  developed 
potential  energy  functions  for  HONO  (Guan  and  Thompson  1989)  and  HCN  (Waite  1984),  with 
modifications  to  accommodate  the  thermochemistry  of  the  system. 

We  have  assigned  the  vibrational  modes  according  to  the  eigenvectors  rather  than  the  values  of  the 
frequencies.  Eigenvectors  of  the  reactant,  transition  state,  and  products  are  in  good  agreement  with  the 
ab  initio  values  (Mowrey  et  al.  1990)  that  are  available,  and  can  all  be  assigned.  Agreement  of  the 
frequencies  for  some  of  the  modes  are  poor.  However,  the  results  of  our  previous  study  (Rice  et  al.  1991) 
showed  that  good  agreement  of  frequencies  at  the  expense  of  the  eigenvectors  can  drastically  affect  the 
reaction  dynamics.  Thus,  we  have  used  the  ab  initio  eigenvectors  for  equilibrium  MN  and  the  five- 
centered  transition  state  as  the  main  fitting  criteria  for  parameterizing  the  PES.  Figure  3  shows  the 
ab  initio  (Mowrey  et  al.  1990;  unpublished)  eigenvectors  (3a)  and  eigenvectors  obtained  from  our  PES 
(3b)  that  correspond  to  the  imaginary  frequency  of  the  transition  state.  The  agreement  is  veiy  good. 
Predicted  frequencies  of  products  by  our  model  PES  are  in  good  agreement  with  the  experimental  values 
except  for  NO  and  NO2,  which  result  from  the  neglect  of  the  modifications  in  the  PES  due  to  the  slight 
changes  in  geometry  and  force  constants  as  these  products  are  formed. 

Mowrey  et  al.  (1990)  estimate  that  the  activation  energy  and  A-factor  for  Reaction  (II)  are 
31  ±  4  kcal/mol  and  10^^  s“^  respectively.  We  refined  the  model  PES  using  this  information  by 
calculating  thermal  rate  coefficients  for  200-1 ,500  K,  and  extracting  Arrhenius  parameters  using  canonical 
variational  transition  state  theory  (Truhlar,  Isaacson,  and  Garrett  1985)  and  the  POLYRATE  5.0.1  set  of 
computer  codes  (Lu  et  al.  1992;  Liu  et  al.  1993).  Adjustments  were  made  to  the  classical  barrier  height 
imtil  the  calculated  activation  energy  was  in  agreement  with  the  Mowrey  et  al.  (1990)  estimate. 

Generalized  transition  state  theory  (GTST)  rate  constants  (Truhlar,  Isaacson,  and  Garrett  1985)  have 
the  form 
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exp(-Vj^£p(s)  /  kjjT), 


(10) 


k®'^(T,s) 


Q^'^CT.s) 
h  Q^^CT) 


where  a  is  a  symmetry  factor,  T  is  the  temperature,  kg  and  h  are  Boltzmann’s  and  Planck’s  constants, 
respectively,  and  s  is  the  distance  along  the  reaction  path  hpom  the  saddle  point.  The  reaction  path  is 
generally  taken  to  be  the  minimum  energy  path  (MEP)  from  the  saddle  point  to  the  reactant  and  product 
geometries.  Q*^^(T,S)  is  the  partition  function  for  the  bound  degrees  of  freedom  at  tiie  generalized 
transition  state  location,  s,  Q*^(T)  is  the  partition  function  of  die  reactant  species,  and  Vj^|£p(s)  is  the  value 
of  the  potential  energy  at  the  generalized  transition  state  location  s.  In  canonical  variational  theory,  the 
rate,  k^^(T),  is  obtained  by  optimizing  equation  (10)  with  respect  to  s.  The  conventional  transition  state 
theory  rate,  k$(T),  corresponds  to  s  =  0. 

We  wish  to  point  out  that  the  rates  calculated  using  equation  (10)  are  not  classical  GTST  rate 
constants,  because  the  partition  functions  for  the  bound  modes  of  the  generalized  transition  states  and  the 
reactant  are  evaluated  quantum  mechanically.  The  reaction  coordinate  motion,  however,  is  treated 
classically.  Quantal  effects  on  reaction  coordinate  motion  are  included  by  using  multiplicative 
semiclassical  ground-state  adiabatic  transmission  coefficients 

kCVT/Y(T)  =  KYkCVT(T)  (11) 

where  Y  is  merely  a  label  that  denotes  the  various  methods  used  to  calculate  these  semiclassical 
transmission  coefficients  (Truhlar,  Isaacson,  and  Garrett  1985).  These  coefficients  generally  account  for 
quantum  mechanical  tunneling.  Because  Mowrey  et  al.  (1990)  did  not  include  tunneling  effects  in  their 
estimate  of  the  activation  energy,  we  used  the  VTST  rates  that  treated  reaction  coordinate  motion 
classically  (k^'^[T])  to  refine  the  PES.  The  interested  reader  is  referred  to  Truhlar,  Isaacson,  and 
Garrett  (1985)  for  a  detailed  description  of  the  theory,  methods,  approximations,  and  additional  references. 

Conventional  and  variational  transition  state  theory  rates  for  Reacticm  (II)  are  shown  in  Table  3. 
Arrhenius  parameters  shown  in  Table  3  were  fitted  over  the  temperature  range  300-1 ,500  K.  The  k^'^(T) 
activation  energy  and  A-factor  are  31.8  kcal/mol  aird  4  x  10^^  s~^  respectively,  in  good  agreement  with 
the  values  estimated  by  Mowrey  et  al.  (1990).  The  rates  that  include  tuimeling  effects  indicate  that 
tunneling  is  not  important  at  temperatures  greater  than  300  K,  but  is  significant  at  lower  temperatures. 
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The  reaction  path  energies  reported  by  Mowrey  et  al.  (unpublished)  were  given  relative  to  the  energy 
at  the  saddle  point  for  (II).  Therefore,  the  reaction  path  energies  were  shifted  accordingly  as  adjustments 
in  the  height  of  the  saddle  point  were  made  during  refinement  of  the  PES  using  VTST  (Truhlar,  Isaacson, 
and  Garrett  1985).  Rgure  4  shows  the  results  of  the  fit  of  the  model  PES  (equation  [1])  to  the  ab  initio 
information.  The  circles  denote  the  points  along  the  ab  initio  reaction  path  (Mowrey  et  al.  unpublished), 
the  crosses  denote  the  model  PES  (equation  [1])  for  the  same  structures,  and  the  solid  line  is  the  reaction 
path  (not  including  zero  point  enei:gy)  of  the  PES  predicted  by  POLYRATE  5.0.1  (Lu  et  al.  1992;  Liu 
et  al.  1993).  The  reaction  path  predicted  by  our  model  is  in  good  agreement  on  the  reactant  side  of  the 
saddle  point  (positive  s).  The  agreement  with  ab  initio  on  the  product  side  (negative  s)  is  not  as  good  as 
for  the  reactant  portion  of  the  path,  but  is  stiU  reasonable.  The  discrepancies  between  the  ab  initio  and 
the  model  PES  energy  curves  are  due  to  differences  in  predicted  product  geometries.  To  illustrate  this, 
the  N-N  (circles),  C-H  (diamonds),  and  H-0  (squares)  bond  distances  along  the  reaction  path  are  shown 
in  Figure  5.  The  lines  in  this  figure  correspond  to  geometries  along  the  model  PES  reaction  path,  and  the 
symbols  correspond  to  the  ab  initio  values  (Mowrey  et  al.  unpublished).  The  model  N-N  and  C-H  bond 
distances  are  in  very  good  agreement  with  the  ab  initio  values  all  along  the  path,  but  the  H-O  distance 
toward  the  products  differs  by  approximately  0.03  A.  This,  and  other  slight  differences  in  product 
geometries  between  the  model  PES  and  ab  initio  calculations  explain  the  energy  differences  for  large 
negative  s  values. 

Figure  6  shows  the  potential  energy  as  a  function  of  the  N-N  distance.  For  each  fixed  value  of  the 
N-N  distance  shown  in  this  figure,  the  remaining  internal  coordinates  have  been  relaxed  to  the  geometry 
at  which  the  first  derivative  of  the  energy  with  respect  to  all  remaining  coordinates  is  zero.  As  expected, 
this  results  in  an  energy  profile  path  for  Reaction  (I)  that  has  no  barrier  beyond  its  endothermicity. 

We  also  foimd  a  saddle  point  leading  to  hydrogen-atom  migration  in  HONO.  The  structure  of  the 
transition  state  species  is  almost  square;  the  HON  angles  are  91.2°,  the  ONO  angle  is  91.0°,  and  the  N-0 
and  H-0  bonds  are  1.2340  and  1.1313  A,  respectively.  The  barrier  to  hydrogen  migration  in  HONO 
predicted  by  this  model  (51.4  kcalAnol)  is  slightly  less  than  the  dissociation  to  HO  +  NO  (52.0  kcalAnol). 
The  actual  migration  barrier  is  likely  considerably  lower  than  our  model  PES  predicts.  We  do  not  expect 
this  higher  migration  barrier  to  affect  the  results  and  conclusions  of  our  planned  dynamics  study,  which 
will  focus  mainly  on  the  primary  decomposition  reacticms  of  CH2NNO2. 
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mblished).  crosses  are  the  energies  predicted  by  equation  (1)  using  the  ab  initio  geometries,  and  the  solid  line  is  the  reaction 
li  predicted  by  equation  (1).  using  POLYRATE  5.0.1  (Lu  et  al.  1992:  Liu  et  al.  1993). 
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Potential  energy  as  a  function  of  N-N  distance  predicted  by  equation  (1).  All  internal  coordinates  of  the  MN  molecule  are  relaxed  to 
the  equilibrium  position  (gradients  are  zero)  for  each  N-N  distance  in  this  figure. 


3.  SUMMARY 


We  have  developed  and  reported  the  analytic  PES  for  MN,  whidi  can  be  used  in  molecular  dynamics 
simulations  of  the  unimolecular  decomposition  of  CH2NNO2.  The  potential  energy  fimction  incorporates 
ab  initio  MCSCF  and  MRCI  calculations  (Mowrey  et  al.  1990;  unpublished)  of  critical  points  for  two 
competing  reaction  paths: 


CH2NNO2  -»  H2CN  +  NO2 


(I) 


and 


CH2NNO2  HCN  +  HONO  .  (II) 

The  bond  dissociation  energy  for  N-N  bond  scission  in  (I)  is  35  kcal/mol;  this  reaction  has  no  back 
reaction  barrier.  The  activation  energy  of  (II)  has  been  estimated  to  be  31  ±  4  kcal/mol  based  on  the 
ab  initio  calculations  (and  corrections)  by  Mowrey  et  al.  (1990).  The  transition  state  for  this  reaction  is 
a  five-centered  cyclic  structure  with  a  classical  barrier  of  37  kcal/mol.  The  barrier  height  corresponding 
to  this  transition  state  was  adjusted  to  predict  a  thermal  activation  energy  of  31.8  kcal/mol,  in  agreement 
with  the  Mowrey  et  al.  (1990)  estimates.  Additionally,  matrices  of  the  energy  second  derivatives  for 
stationary  points  and  points  along  the  reaction  paths  that  were  calculated  by  Mowrey  et  al.  (unpublished) 
were  used  in  parameterizing  the  PES. 

Thermal  rates  for  (II)  were  calculated  using  Variational  Transition  State  Theory  (Truhlar,  Isaacson, 
and  Garrett  1985)  and  used  to  adjust  the  PES  until  they  were  in  agreement  with  the  Mowrey  et  al. 
estimates  (Mowrey  et  al.  1990).  Tunneling  effects  in  (II)  are  significant  only  at  low  temperatures,  and 
do  not  contribute  appreciably  to  the  rate  of  decomposition  at  the  energies  considered  in  this  study. 
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